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ABSTRACT 

We study the outer-shock structure of the oxygen-rich supernova remnant G292.0-I-1.8, usmg a 
deep observation with the Chandra X-ray Observatory. We measure radial variations of the electron 
temperature and emission measure that we identify as the outer shock propagating into a medium 
with a radially decreasing density profile. The inferred ambient density structure is consistent with 
models for the circumstellar wind of a massive progenitor star rather than for a uniform interstellar 
medium. The estimated wind density (njj = 0.1 ~ 0.3 cm~^) at the current outer radius (~ 7.7 pc) 
of the remnant is consistent with a slow wind from a red supergiant (RSG) star. The total mass of 
the wind is estimated to be ^ 15 — 40 Mq (depending on the estimated density range), assuming that 
the wind extended down to near the surface of the progenitor. The overall kinematics of G292.0-I-1.8 
are consistent with the remnant expanding through the RSG wind. 

Subject headings: ISM: individual objects (G292.0-I-1.8) — ISM: supernova remnants — shock waves 
— stars: mass-loss — X-rays: ISM 



1. INTRODUCTION 

Studying the environments in which supernovae (SNe) 
explode and then subsequently evolve is crucial to under- 
stand how SNe feed energy and matter into the Galactic 
ecology. This is particularly important for core-collapse 
SNe which are believed to explode in a complicated envi- 
ronment caused by the late stage evolution of a massive 
progenitor star. The study of young core-collapse super- 
nova remnants (SNRs) interacting with their circumstel- 
lar media is important in this regard. The evolution of 
the SNR is not only governed by the medium into which 
it is expanding but also by the details of its explosion. 
Therefore, the study of young core-collapse SNRs may 
also provide valuable information about their explosion. 

The environment of a massive star at the time of 
the explosion is complicated by its evolutionary history. 
Stars of initial mass < SSM© are likely to explode as 
red super giants (RSGs), while more massive stars end 
as Wolf- Rayet (WR) s tars with a possible earlier RSG 
phase (|Chevaliei]l2005l ). For a simplistic case where the 
stellar wind had a constant velocity and a constant mass 
loss rate, the circumstellar density profile of the freely 
flowing unshocked wind is given by p oc r^^ where p is 
the density of the wind at the distance r from the star. 
The structure and evolution of a SNR expanding in such 
a medium differs from t hat of a SNR e xpanding in a 
uniform density medium. iChevalieil ()1982j ) derived a self- 
similar solution for the 1-dimensional spherical expansion 
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of SN ejecta whose density profile is given as a power-law 
{p cx r^") in a medium with a radial density dependence 
{p cx r^^). The evolution beyond the RSG wind is more 
complicated (i.e., similarity solutions fail to apply) and 
numerical studie s are often require d for a comprehensive 
description (e.g.. lDwarkadasll2005l ). A Wolf-Rayet phase 
expected in more massive s tars further c omplicates its 
envir onment (e.g., Chevalier fc Lian"^ 1198 9: Dwarkadai 
[200l . 

The circumstellar interaction of SNe has been ex- 
tensively studied, especially in the radio band (see 
iWeiler et al.l I200"2L for a review). Such an interaction 
plays a criti cal role in interpreting observations of 7-ray 
bursts fe.g.. iCampana et al.ll2006f) . On the other hand, 
observational studies of core-collapse SNRs expanding 
within a circumstellar wind have not been well estab- 
lished so far. Among more than 260 Galactic SNRs, only 
a few known core-collapse SNRs sho w circumstantia l evi- 
dence of the shock-wind interaction (lChevalieill2005L and 
references therein) , among which the most notable exam- 
ple is Gas A. The overall characteristics of Gas A are well 
described by a model of a supernova interacting with a 
RSG wind with a mass loss rate of ~ 2 x 1O~^M0 yr~^ for 
a wind velocity of 10 kms~^ (|Chevalier fc Oishill2003[ ). 

Another well-known example of a young core-collapse 
SNR is G292.0+1.8. G292.0-H. 8 is one of the 
three known "oxyge n-rich" SNRs (jGoss et al.l 119791: 
iMurdin fc Clarklll979D in the Galaxy, along with Gas A 
and Puppis A. These SNRs show optical emission from 
fast-moving {v > 1000 km s~^) 0-rich ejecta knots, which 
are suggested to be synthesized in th e core of a massive 
star (> 10 Mq; e.g. iBlair et"alll2000[) . The existence of 
the central pulsar (P S R J1124-5916 ) and i t s wind nebula 
llHughes e t all 120011: iCamilo et all 120021: iHughes et all 
2003; Gae nsler fc Wallace! I2003| ) clearlv demonstrates 
that G292.0-I-1.8 originated from a core-collapse SN of 
a massive star. A progenitor mass of ~ 20 — 40 Mq 
is suggested from the observed metal abundance struc- 
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ture of an X-ray emitting ejecta (iHughes &: Singhlll994l : 
IGonzalez fc Safi-Harbl \2oM IPark et alJl2004D . Previous 
X-ray studies showed an X-ray bright central belt-hke 
structure, which hkely c orresponds to shoc ked circum- 
stellar medium (GSM) ()Park et all I2OO20 . iPark et alJ 
()2002[ ) suggested that the belt traces the enhanced 
mass loss in the progenitor star's equatorial plane. 
G292.0-I-1.8 also shows faint X-ray emission along the 
outer boundary of the remnant, which corresponds to 
a sharp-edged perimeter seen in the radio morphology 
(jGacnslcr & Wallace 2003). This faint X -ray emission is 
produced by the blast wave propagating into the ambient 
medium. It shows a simple, nearly circular morphology 
all around the SNR without showing significant substruc- 
tures and is mostly free from shocked metal-rich ejecta 
emission. Thus, this outer region of G292.0-(-1.8 is use- 
ful for the detailed study of the GSM and its interaction 
with the SN blast wave. 

In this paper, we present the results from our analy- 
sis of shocked GSM in G292.0-M.8 using our deep ^500 
ks Gh andra observation first presented by iPark et all 
(|2007| ). which provides an excellent opportunity for the 
spectral study of this faint emission from the outer 
boundary of the SNR. We describe the Ghandra obser- 
vation and the data reduction in § [5] In §§ [3] and SI 
we investigate radial and azimuthal variations of the X- 
ray spectral properties of shocked GSM to show that the 
remnant is expanding in a circumstellar wind. The impli- 
cations of our results on the nature of the circumstellar 
wind and the kinematics of the remnant are discussed in 
§[5] §[5] summarizes our results. 

2. OBSERVATIONS AND DATA REDUCTION 

We observe d G292.0-H.8 on 2006 September 13 - Oc- 
tober 20 (Park et al.l I2007D using the Advanced GGD 
Imaging Spectrometer (AGIS) on board Chandra. The 
entire remnant (~ 9' in diameter) is covered with AGIS-I 
(17' X 17' field of view) and the pointing was selected to 
place the pulsar (PSR J1124-5916) near the aimpoint. 
The total effective exposure time is 506 ks after the 
data reduction, with the exposure of individual obser- 
vations ranging between ~ 40 and 160 ks (Table [ij. The 
level 1 event files were reprocessed to create new level 
2 event files. We applied parameters of the standard 
Ghandra pipeline process, except that we applied the 
VFAINT mode background cleaning and turned off the 
pixel randomization. GIAO 3.4 and GALDB 3.2.3 were 
used for all the reprocessing and analysis. We examined 
the overall background light curve for periods of high 
background. Only minor increases above the quiescent 
background level by a factor of < 2 were noticed and 
their total duration is less than 1% of the total exposure. 
Therefore, we did not apply any light curve filtering. 

We checked the relative astrometry between individ- 
ual observations using the position of the pulsar (PSR 
J1124-5916). We found less than 0.1" deviations which 
are within typical on-axis astrometric uncertainties of the 
AGIS. For the purpose of the imaging analysis, event files 
were reprojected to a common tangent point (that of Ob- 
sID 6679). We generated exposure maps for each ObsID 
with the dead area correctiorQ applied. For the spectral 
analysis, spectrum and response files were generated for 



TABLE 1 

Summary of Chandra observations 



ObsID 


Date 


Exposure 


RoUAnglct 






[ks] 


[°] 


6680* 


9/13/2006 


39 


180 


6678 


10/2/2006 


44 


157 


6679 


10/3/2006 


153 


157 


8447 


10/7/2006 


47 


157 


6677 


10/16/2006 


159 


140 


8221 


10/20/2006 


64 


140 



t Roll angle describes the orientation of the 
Chandra instruments on the sky. 
t ObsID 6680 was conducted with 2 ACIS-S chips 
in addition to 4 ACIS-I chips, and affected by 
telemetry saturation, reducing the net exposure 
time by ~ 25%. All later observations used only 
ACIS-I chips, and thus were not affected by the 
telemetry saturation. 

individual ObsIDs and then merged, weighting by the ex- 
posure time. Out-of-time event^ from the bright SNR 
emission along the GGD readout directions are present, 
which in principle could affect our spectral analysis. For 
the southeastern region, where the SNR emission is faint, 
the estimated out-of-time event counts are 7% of the to- 
tal counts (source and background). However, the out- 
of-time events contribute uniformly along the GGD read- 
out direction and we find the local background subtrac- 
tion removes more than 90% of out-of-time events in the 
source spectrum. The effects of out-of-time events to the 
background-subtracted spectrum are thus negligible. 

3. RADIAL PROFILE OF THE OUTER SHELL OF 
G292.0+1.8 

3.1. Spectral Fitting 

The X-ray images of the deep Ghandra observation 
clearly reveal the faint outer blast wave all around the 
remnant (Fig. [1} . While the northern and western outer 
boundaries are well-defined and relatively bright, that of 
the southeast is faint. Although a slight limb-brightening 
appears to be present in some parts of the northwestern 
boundary, in general, the surface brightness of the outer 
boundary does not show a significantly limb brightened 
morphology. The radial brightness variation is rather fiat 
and even increasing inward in some regions (e.g., south- 
east). Given the core-collapse origin of the remnant, we 
consider this overall morphology in the framework of the 
SNR expanding in its circumstellar wind. 

For a quantitative analysis, we investigate the radial 
variation of the X-ray spectral characteristics near the 
outer boundary. This approach is motivated by the fact 
that the gas behind the forward shock preser ves the in- 
forma tion of the medium it has swept up. iGhevalierl 
(|1982f) found self-similar solutions for the interaction of 
expanding stellar ejecta {p oc r^", n > 5) with an ex- 
ternal medium (p cx r~*): a remnant expanding in a 
uniform ambient medium (s = 0) would result in an 
inward-decreasing density and an inward-increasing tem- 
perature profile for the shocked ambient gas, while a rem- 
nant expanding in a medium with a wind profile (s = 2) 
would show an inward-increasing density and an inward- 
decreasing temperature. Therefore, by studying the cur- 
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rent radial profile of the density and temperature, we 
can infer whether the SNR is currently expanding into a 
density gradient. 

We selected two representative regions in the north- 
west (NW) and southeast (SE). We examined various 
narrow band images and X-ray color images of the rem- 
nant and selected regions where the contamination by 
ejecta emission is negligible and the area is sufficiently 
large to study the radial spectral variation. Each re- 
gion (NW and SE) is divided into several radial subre- 
gions to track the radial spectral variation (Fig.[T]). The 
photon statistics for each subregion are ^ 3000 — 5000 
counts in the 0.5-8 keV band. The extracted spectrum 
for each subregion is fitted with a single nonequilibrium 
ionization ( NEI) plane-parallel sho ck model (vpshock in 
Xspec vl2. iBorkowski et al.ll2001| ). The fits with vari- 
able metal abundances (O, Ne, Mg, Si, S and Fe were 
fitted while the other species w ere fixed at their solar 
values (| Anders fc Grevessill989[ )) are statistically better 
{x^ 1^ ^ 1-0) than those with metal abundances frozen 
{x^/v ^ 1.2) at solar values. The fitted metal abun- 
dances are generally subsolar, confirming that there is 
no significant ejecta contamination in the selected re- 
gions. The derived spectral parameters from the indi- 
vidual subregional fitting are not well constrained due to 
the large number of free parameters involved in the fit. 
As the abundance pattern in each radial region is similar 
within statistical uncertainties, we simultaneously fit all 
subregions in each of the NW and SE regions with their 
metal abundances tied. The variation of the absorbing 
column density when fitted independently is statistically 
insignificant, and thus we tie the absorbing column den- 
sities among subregions, as well. We note that our fits 
result in a relatively higher iVn of ~ 5 — 6 x lO^-'^cm^ than 
TVh ~ 3.2 X lO^^cm'^ from the spectral analysis of the 
pulsar (|Hughes et al.ll2001[ ). Fixing TVr at 3.2 x lO^^cm^ 
increases the temperature (by ~ 20%) and decreases the 
normalization (by ~ 30%) in a systemic way and does not 
fundamentally affect our results. The above approach 
provides acceptable fits (example spectra and the best- 
fit models are shown in Fig. [2]) for individual subregions 
in each of the SE and NW regions with well constrained 
spectral parameters. The fitted parameters are consis- 
tent with those from independent fits within statistical 
errors. We utilize these fit parameters for the following 
radial analysis. The fit results for regions NW and SE 
are summarized in Tables [2] and |3l 

3.2. Model Comparison 

In Fig.[3l we plot the spectral variation as a function of 
the radial distance of each subregion from the geometric 
center of the remnant as determined by the radio con- 
tinuum image (a2nnn, <^2nn n = ll'' 24™ 34!8, -59°15'52'.'9, 
iGaensler fc Wallacg|2003[ ). where the radial distance of 
each subregion is normalized by the distance to the for- 
ward shock in that region. As these regions are located 
in the outer parts of the remnant, the uncertainty of the 
remnant center does not make a significant effect on the 
results presented in this section. The radial variations 
of the spectral parameters are consistent between the 
NW and SE regions (Fig. [S]) . The emission measur^^ in- 

^ We use the conventional emission measure (EM = / rienn dl) 
where the fitted volume emission measure (J nenjj dV) is normal- 



counts cm ^ s ^ pixel ^ 
0.5 X 10-' 1 X 10"' 




Fig. 1.— The 0.3 - 10 keV band Chandra ACIS image of 
G292.0-I-1.8 with regions used in the radial spectral analysis over- 
laid. The size of each individual region is adjusted to have a similar 
number of photon counts while avoiding any local features. The 
image is exposure-corrected and has a pixel size of 0'.'492. The grey 
scale is saturated in the bright filamentary structures near the cen- 
tral regions of the SNR to emphasize the faint outer shock. T he 
radio geometrical center of the SNR IGaensler fc Wallace! [2003) is 
marked as a white X. 
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Fig. 2. — The Chandra ACIS spectra of the SE sector. Spectra at 
different radial distances (SEl, SE3, SE5) are shown. The spectra 
of SE3 and SE5 are arbitrarily scaled for the purpose of display. 
The solid lines are the best-fit models. 

creases inward, while the electron temperature does not. 
The ionization timescale {net) shows a similar trend as 
the emission measure. 

We overlay model predicted radial distributions of the 
temperature and e mission measure in Fig. [31 where the 
radial profiles from lChevalieii (Il982f) are pr ojected assum- 
ing a spherical geometrv. iChevaliCT ()1982[ ) provides mod- 
els of density, pressure and velocity of the g func- 
tion of the radial distance from the center. The model 
emission measure in Fig. [3] is calculated by integrating 
the density square along the line of sight. The tem per- 
ature is not explicitly provided in iChevalieil (|1982D , and 
we have calculated it from the thermal pressure and gas 
density. The projected temperature in Fig. [3] is taken as 
an average temperature along the line of sight weighted 

ized by the region area. 
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by density squared. While we ignore the variation in 
the intrinsic emissivity of the gas at different temper- 
atures, the variation is not significant (< 30 %) in the 
temperature range of 0.6 - 0.9 keV. Models with four dif- 
ferent sets of n and s are plotted. The solutions of s = 2 
show distinguished radial variations from the solutions of 
s — 0, which reflect their different radial trends in den- 
sity and temperature as briefly described in § 13.11 The 
projected temperature of the shocked ambient medium 
of s = 2 solutions shows a decrease toward the contact 
discontinuity (CD), in contrast to the increase in s = 
solutions. The emission measure of s = 2 solutions mono- 
tonically increases toward the CD. The emission measure 
of s = solutions initially increases due to the geomet- 
rical effect, and then starts to decrease. These overall 
trends of the emission measure and temperature profiles 
are dominated by the ambient density structure (s = 
vs. 2) rather than by the ejecta profile. The primary 
effect by the power law index n of the ejecta is on the 
distance between the forward shock front and the CD: 
e.g., a steeper ejecta profile makes the CD close to the 
forward shock. Fig. [3] demonstrates that the general ten- 
dency of the observed temperature and emission measure 
are well described by models with s = 2 (the solid and 
dashed lines). The model of s = 2, n = 7, in particular, is 
in good agreement with the observations, where our fits 
show reduced of 0.5 ^ 1 in both SE and NW regions. 
On the other hand, s — models (the gray and the dot- 
ted lines) show significant discrepancies especially in the 
SE region with a reduced of ~ 3. The observed radial 
structures of the temperature and the emission measure 
are inconsistent with the remnant expanding in a uniform 
medium where the postshock inner region is expected to 
result in a lower density and a higher temperatu re than 
the r egion immediately behind the blast wave (Chevalieii 
[198^ . 

We note that the temperature from the X-ray spec- 
tral fit is the electron temperature, while the temper- 
ature from the model is the ion temperature. These 
two can differ considerably in non-radiative shocks faster 
than 1000 kms~^ (e.g.. iGhavamian et aIll2007L and ref- 
erences therein). Assuming the SNR age o f 3000 yr 
(jChavamian et all 120051 : iWinkler et all 120091 this age 
is adopted throughout this paper), the average blast 

wave velocity would be 2500 km s~^ for a remnant ra- 
dius of ~ 7.7 pc (see § 15. 2p . and the blast wave should 
have been faster in the past. The observed ratio be- 
tween the electron temperature (Te) and the proton tem- 
perature (Tp) immediately be hind t he shock is known 
to vary with shock velocities (jGhav amian et al. 2007). 
But the detailed variation of Te/Tp in the shocks of 

Vs ^ 2500 km s~^ is unclear and we assume a constant 
value of To/Tp across the radius of the postshock re- 
gions used in our analysis (the actual value does not 
matter for our analysis). We note that Coulomb equi- 
libration in the downstream region will eventually equi- 
librate ions and electrons, and the electron temperature 
will increase. This would increase the downstream elec- 
tron temperature. Qualitatively, the modeled radial pro- 
file of the electron temperature would thus be fiatter (in 
case of the s = 2 models) than that of the ion tempera- 
ture, which would still b e consistent with our data. 
The analytic model of iChevalieri (|1982| ) is a similarity 



solution assuming a power law distribution of the ambi- 
ent density and ejecta, which may not be entirely ade- 
quate for actual SNRs. For example, the s = 2 models 
have the initial wind density profile of infinite at r = 
and predict an infinite density and zero temperature at 
the CD. However, we believe that this solution is a good 
approximation to the real case for the regions away from 
the CD. While the density distribution of the ejected 
material depends on the structure of the progenitor star, 
a power law envelope surrounding a flat (or shallower) 
core i s a reasonable approxim ation (Chevalier & Sokea 
119891 : iMatzner fc McKeel [l999[ ) . The assumption of the 
power law density of ejecta will eventually break down 
as the reverse shock reaches the core. While this may 
significantly modify the structure of the shocked ejecta 
material, the effect on the structure of the shock e d am - 
bient gas would not be significant. The lChevaliei] ()1982[ ) 
models do not include the effects of cosmic ray acceler- 
ation in the shock. While the efficient cosmic ray ac- 
celeration may modify the underlying shock structure, 
the overall tre nd of the radial variation is likely to be 
stiU vahd (e.g.. iDecourchelle et al.ll2000l ). Therefore, we 
conclude that the observed trend of the decreasing tem- 
perature and increasing emission measure of the shocked 
ambient gas toward the CD is a robust indicator that 
the blast wave is currently expanding into a radially- 
decreasing density profile (s = 2) rather than a constant 
density (s = 0). This result provides direct observational 
evidence that the blast wave of G292.0-1-1.8 is currently 
propagating through the stellar winds produced by its 
massive progenitor. On the other hand, our results may 
not constrain the power law index of the ejecta (n). The 
apparent effects of n is the change in the relative dis- 
tance from the CD to the outer boundary. However, this 
distance can be varied by other physical processes which 
are not properly included in the model, e.g., cosmic ray 
accelerations. Also, the propagation of the reverse shock 
into the shallow ejecta core may result in the decrease 
of the effective power-law index. Therefore, while the 
observation shows better agreement with the model of 
n = 7 than the model of n = 12, this may not be in- 
terpreted in favor of the n = 7 model over the n = 12 
model. 

4. AZIMUTHAL VARIATION 

Fig. 21 shows regions we chose for the azimuthal study. 
We selected narrow regions (thickness 10" ^ 20") along 
the outer boundary of the remnant. We avoided regions 
in which the contamination from the ejecta emission is 
significant We extracted the spectrum from each region 
and fitted it independently with the vpshock model. The 
photon statistics in these regions are a few thousand 
counts which are similar to those in subregions used in 
the radial analysis. The fitted spectral parameters are 
poorly constrained when the metal abundances are freely 
varied, and we fix the metal abundances to the average 
values as obtained from the radial analysis (NW and SE 
regions). The fitted absorbing column density varies be- 
tween 5.5 and 7.0 x 10^^ cm~^. However, the Nh vari- 
ation is statistically insignificant, e.g., the fits do not 
significantly change when Nh is fixed at the mean value. 
Thus, in the following discussion, we use results from the 
fits where we fix Nh at their mean of 6.2 x 10^^ cm~^. 

The azimuthal variations of the electron temperature 
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TABLE 2 

Fitted A^h and metal abundances for the regions used for the radial analysis. 



Sector 




O 


Ne 


Mg 


Si 


S 


Fe 


NW 
SE 


c: O + 0.4 

^■^-o.s 


67+°-^^ 
'-'•"'-o.io 

36+° " 


n 75+0-0'' 

U-'0_o.07 


n .:>r+0.09 


.0+0.03 
"•^'-'-0.03 


U.OO_Q 24 
12+0" 


24+0 06 
07+0 03 

^■^'-0.02 



Note. — These parameters are fitted after being tied among subregions. Results for spectral parame- 
ters v aried freely are presented in Table[3] Metal abundances are with respect to solar l|Anders Hi Grevessd 
[19891 ). Throughout the paper, errors are shown in 90% confidence range unless explicitly specified other- 
wise. 



TABLE 3 

Spectral parameters for NW and SE regions. 



Region Distance 




log net 


E.M. 


['] 


[kcV] 


[cm~3s] 


[1016 cm-5] 



NWl 


4.33 


0.86 


NW2 


4.21 


0.79 


NW3 


4.11 


0.79 


NW4 


4.03 


0.56 


NWS 


3.94 


0.71 


NW6 


3.85 


0.64 


SEl 


4.30 


0.88 


SE2 


4.08 


0.83 


SE3 


3.87 


0.70 


SE4 


3.65 


0.73 


SE5 


3.46 


0.66 



-1-0.11 
-0.09 
-1-0.08 
-0.06 
+0.11 
-0.03 
+0.12 
-0.03 
+0.15 
-0.06 
+0.11 
-0.04 



+0.16 
-0.10 
+0.12 
-0.07 
+0.11 
-0.11 
+0.10 
-0.04 
+0.10 
-0.09 



10.41 
10.74 
10.75^ 



+0.10 
0.06 
+0.14 
-0.13 
-+0.10 
-0.13 
11 34+0.18 

10 Q3+0-30 
xu.»o_Q 20 

11 05+0-13 



10.69t°:ig 

11 13+0-15 
ll-lO-O.lO 
11 51+0-0'? 

07+0.17 

+0.21 
0.18 



37+0.10 
'J-'" -0.07 
1+0.15 
-0.11 

0.68+°" 



0.64-* 



11.32 



1.01 
0.75 
1.37 

0.10 
0.16 
0.23 
0.30 
0.36 



0.15 
+0.36 
0.36 
+0.34 
0.27 
+0.36 
0.36 



+0.03 
0.02 
+0.05 
0.04 
+0.04 
0.04 
+0.08 
0.07 
+0.12 
0.15 




Note. — Regions are listed on the decreasing order of their 
distance from the center; i.e., the top row in each region repre- 
sents the outermost region. Based on the simultaneous fit for 
each region, x^/i' = 3039/3042 and 2415/2534 for NW and SE 
regions, respectively. 



TABLE 4 

X-RAY spectral PARAMETERS OF THE REGIONS USED IN 
THE AZIMUTH AL ANALYSIS. 



P.A.t 








nil 


1 


og net 






n 




[keV] 




cm~3] 


[lOlO cm-3s] 




13 





Sfi+O-OT 
°tl-0.07 


1 


0(5+0.08 
"•3-0.06 


10 


00+0.15 
■J"-0.12 


1 


06 


59 





76+0-" 
'"-0.22 





qr,+0.27 
^•^-0.08 


10 


7Q+0.61 
'^-0.35 





89 


126 





7-1+0.09 
'1-0.10 
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t Position angles to each region are measured east of north 
from the geometrical center. 

t For all individual spectra, the degree of freedom {u) is 
372. 

and the hydrogen number density are presented in Ta- 
ble m and Fig. [S] The densities are estimated from the 
volume emission measures. The volume corresponding 
to each region is calculated by pixel-by-pixel integra- 



FlG. 3. — Observed electron temperatures and emission measures, 
plotted as a function of the radial distance from the SNR cen- 
ter, for NW and SE regions. The radial distances are normalized 
by distance from the geometric center of the SNR to the blast 
wave (Rhw) ™ each direction. Lines are model predictions from 
IChevalieil J1982I : see § |3] for more details) for models of different 
set of parameters, where s and n are power-law index for the den- 
sity profile of the ambient medium and of the expanding ejecta, 
respectively. Cusps in these lines mark the location of the con- 
tact discontinuity in the model. The temperatures and emission 
measures of the model predictions have been scaled to match the 
observed values. 

tion of the path-lengths over the projected area. In this 
calculation, we assumed a sphere centered at the geo- 
metrical center of the rem nant at a distance of 6 kpc 
([Gaensler fc Wallacil2003[ ). The radius of the sphere is 
separately determined for each region as the projected 
distance of the farthest (from the geometrical center of 
the remnant) pixel in the region. The uncertainties on 
the electron temperature are too large to constrain any 
possible variations around the SNR boundary. On the 
other hand, in the SE regions (PA - 120° - 150°), the 
density shows a significant drop by a factor 2 from the 
average (~ lcm~^). The gas density shown in Fig. [5] 
is the postshock hydrogen number density. Assuming a 
compression factor of 4 from a strong gas shock, this cor- 
responds to a preshock gas density of ^ 0.1 cm"'^ toward 
the SE and 0.2 — 0.3 cm^"^ in other regions. 

5. DISCUSSION 

Our observation provides direct evidence of 
G292.0-1-1.8 interacting with its GSM wind, which 
is consistent with the core-collapse origin of the remnant 
(§1). Our results also can help reveal the nature of the 
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Fig. 4.— The 0.7 - 0.84 keV band Chandra ACIS image of 
G292.0+1.8. The energy band has been selected to emphasize the 
equatorial belt. The regions along the remnant's outer boundary 
used for the azimuthal analysis are marked as black boxes with the 
corresponding position angles (north through cast) indicated. The 
radio geometrical center of the SNR tGaensler fc Wallacei2003il is 
marked as a white X. The position of the pulsar J1124-5916 is 
indicated by a white arrow. 
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Fig. 5. — The azimuthal variation of the electron temperature 
and the postshock hydrogen density. Position angles are measured 
from the geometrical center as counter-clockwise from the north. 
The error bars in the position angle represent the angular extent 
of each region in Fig. |4] 

progenitor star by the observed wind properties. In the 
foUowings, we discuss the nature of the wind with which 
G292.0-I-1.8 is currently interacting and the overah 
kinematics of the remnant. We also briefly discuss the 
possible evolutionary paths of the progenitor star. 

5.1. The Progenitor's Stellar Wind 

The expansion of SN blast waves within circumstel- 
lar winds has been extensively studied with SNe obser- 
vations. The shocked CSM emits in radio continuum 
and also in X-rays. The mass loss rates of the progeni- 



tor stars, estimated by modeling their light curve, are ~ 
10"^ Mq yr-i for type Ib/c and 10""^ - lO"*^ Mq yi'^ 
for type II SNe (ISramek fc Welled [200l . The circum- 
stellar interaction of SNRs, on the other hand, has not 
been well established, although there has been some cir- 
cumstantial evidence (iChc valicr 2005,) . 

The radial profiles of the postshock temperature and 
density in G292.0-I-1.8 argue that the remnant is expand- 
ing within the circumstellar wind. A slight limb bright- 
ening is seen in a few regions in NW, and we attribute it 
to a local effect, e.g., the shock is probably encountering 
locally dense regions in the wind. The overall dynamics 
of the remnant would not be significantly altered by such 
a small local effect, and we consider that the SNR shock 
has not yet reached the termination shock of the wind. 
The radial profile of the wind density with a constant 
mass loss rate of M and a wind velocity of Vw is given 
by p = M /Anr^Vw. Adopting the geometrical center of 
the remnant as the wind center, where we estimate an 
average shock radius of rf, = 7.7 dg pc, we estimate 



M = 2xl0-M 

VlOkms"^ 



f_nB__^ 
vO.l cm^3y 



v7.7pc 



Mq yr- 



-1 



A slow wind velocity (w^ ^ lOkms ^) appropriate 
for the RSG winds implies a mass loss rate M = 
2 — 5 X 10~^MQyr^^ for the observed density range 
(riH — 0.1 — 0.3 cm~'^), which is consistent with the ob- 
servations of type II SNe. On the other hand, assuming 
a fast {v^ > 1000 km s"""^) Wolf-Rayet wind results in 
M > lO~'^M0yr~^, a few orders of magnitudes larger 
than what is inferred from type Ib/c SNe observations. 

The integrated mass of the wind swept up by the rem- 
nant, assuming the wind extends down to r = 0, is 



M,„ 15 



0.1 cm 



n 



7.7pc 



Mr. 



Accounting for the azimuthal density variation in the 
postshock regions, we estimate the mass range of Mw ~ 
15 — 40 M0. There are further uncertainties involved in 
this estimate due to the various assumptions we made. 
Our density measurement scales as inverse square root of 
the assumed distance to the remnant (p cx d^^/^), thus 
Mu, cx d^/^. For example, a slightly smaller distance of 5 
kpc (which is within unc ertainties in the distance mea- 
surements to G292.0-H.8. IGaensler fc Wallacel[200l will 
reduce the mass estimate by ~ 40%. Another source of 
uncertainties could be our simple assumption on the rem- 
nant geometry, a sphere centered at the geometric center 
of the remnant. If the wind is clumped, a local density 
fluctuation may have led us to overestimate the average 
density, as the observed emission measure is likely biased 
by the high density regions. The mass estimate will also 
be reduced if the shock compression ratio is higher than 
4 due to efficient CR acceleration, while this effect may 
not be significant as there is no evidence for nonther- 
mal X-ray emission in G292.0+1.8. We assume that the 
mass loss remained constant until the moment of the ex- 
plosion. However, observations of some SNe suggest that 
their mass loss rates have decreased with time before the 
explosion (e.g.. Ivan Dvk et 311119941) . Also, the progeni- 
tor star might have evolved through different stages (such 
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as Blue Super Giant or Wolf-Rayet star). 

Theoretical studies of the stellar evolution suggests 
that stars of less massive than 35 Mq end their life dur- 
ing the RSG phase, while more massive stars evolve to 
WR phase with a possible earlier RSG phase. Given 
that the wind is likely a dense RSG wind, the progenitor 
star of G292.0-I-1.8 might have a initial mass less massive 
than 35 Mq. In this regards, the wind mass as large as 
40 Mq may not be realistic for stellar evolutionary mod- 
els. However, we note that there could be uncertainties 
in the modeling of the late stage of the stellar evolution, 
which is very sensitive to mass loss rate of the star, as 
the stellar mass loss rate is rather poorly constrained 
by observations. On the other hand, previous nucle- 
osynthesis studies of G292.0-I-1.8 using X-ray emission 
from metal -rich ejecta suggested a progenitor mass of 
20-4 M^ (IHug hes fc Singhll994l :IGonzalez fc Safi-"Ha7bl 
[200llPark eTaL 20041. while these progenitor mass es- 
timates were limited by nondetectio n of the explosiv e 
nucleosynthesis products in this SNR (jPark et al.ll2004[ ). 
Considering that there are large uncertainties on both 
estimates of the wind and the progenitor masses, our es- 
timate of the wind mass is in plausible agreement with 
a general picture: the progenitor of G292.0+1.8 was a 
massive star which lost a significant amount of its mass 
during the pre-SN evolution. 

Our azimuthal analysis in § |4] suggests an azimuthal 
density variation. Since the outer boundary of the rem- 
nant has unlikely reached the termination shock of the 
RSG wind, the observed density variation would repre- 
sent the variation in the wind density. The wind density 
in the SE boundary of the remnant is smaller than those 
in other regions by a factor of ^ 2. While this could 
be due to some local effects, the wind might have been 
asymmetric. Alternatively, the center of the wind might 
have not been at the geometric center of the SNR. The 
asymmetric SN explosion might have caused the south- 
eastern part of the remnant to expand more rapidly than 
in NW, and the azimuthal density variation could be due 
to an offset of the progenitor s tar from the SNR's geo- 
metric center. Recently, iWinkler et al.l (2009) studied 
proper motion of [O III] ejecta and located the estimated 
expansion center of ejecta knots, which turned out to be 
consistent with the SNR's geometric center. However, 
we note that the observed proper motions of many indi- 
vidual optical knots show significant deviations from the 
best-fit model which indicated the explosion center at 
the geometric center of the SNR. Thus, the origin of the 
azimuthal density variation along the outer boundary of 
G292.0-I-1.8 is not conclusive. 

5.2. Supernova Remnant Kinematics 

The morphology of the remnant, especially when it 
is young, traces its explosion characteristics. The kine- 
matic properties of G292.0-I-1.8, such as the age and 
the explosion energy, have been estimated in the pre- 
vious works. However, G292.0-I-1.8 was previously as- 
sumed to be expanding in a un iform medium instead 
of the circumstellar wind (e.g. , IGonzalez fc Safi-Harbl 
I200I IGaensler fc Wallacg l2003l) . although G292.0-hl.8 
is believed to be a core-collapse SNR. In fact, detailed 
observational studies of young core-collapse SNRs ac- 
counting for the evolution in a circumstellar wind ma- 
terial have been conducted for only a limited number 



of remnants such as Gas A (jGhevalier fc Oishl 120031 : 
ILaming fc Hwan^ [2003). Here, we estimate kinematic 
properties of G292.0-fl.8, explicitly accounting for its 
evolution in the circumstellar wind. 

For simplicity, we assume a spherically symmetric 
explosion. Like the SNR evolution in a homoge- 
neous medium, the initial evolution of the remnant 
is dominated by the characteristics of the ejecta (the 
ejecta-dominated phase) and makes a transition to the 
Sedov- Taylor (ST) phase as the s wept-up mass be- 
comes larger than the ejecta mass. i Truelove fc McKo^ 
(|1999l TM99 hereafter) discussed a unified solution 
which describes a spherically symmetric interaction of 
the ejected stellar material, which has a flat core sur- 
rounded by a power-law envelope, with an ambient 
medium. TM99 f ocused on the evol u tion i n the homo- 
geneous medium. ILaming fc Hwan^ (2003, LH03 here- 
after) extended TM99's work to evolution in a stellar 
wind. The characteristic age (to) and size (xo) of the 
remnant, which mark the approximate transition from 
the ejecta-dominated (ED) phase to the ST phase, are 
given by 



to = 5633Mf/^£:5Y^^ 



yr 



xo = 40.74Mej(nHi?fc')"' 



pc 



for the evolution of the remnant in the pre-supernova 
wind (from Eq. A3 and A4 of LH03); nn is a hydro- 
gen number density of the wind at the radius of the 
blast wave (Rb), -E51 = ergs, E is the explo- 

sion energy and M^j is the ejected mass from the su- 
pernova, and Mej, R^ and nn are in units of Mq, pc 
and cm~^, respectively. If we adopt an average radius 
of Rb — 7.7 pc from the geometrical center and the low 
density of nn = 0.1 cm~'^ that seems to be preferred by 

^ /'2 1/2 

stellar evolution model, we obtain to = ^b^dM'^'j^ E^^^' ^ 
yr and xq = 6.9Mej pc for G292.0-I-1.8. Assuming M^j 
of several Mq as in Gas A and for a canonical explosion 
energy of i?5i = 1, we estimate to ^ 10^ yr, which is 
larger than the remnant age (3000 yr) (technically, the 
tra nsition time to the ST pha se (tsr) can be smaller than 
to, T ruelove fc McKed [T999f) . Thus G292.0-H1.8 might 
still be in its ED phase or in its transition to the ST 
phase. And the ED phase evolution would be better fit 
to describe the current kinematics of the G292.0-I-1.8. 
The evolution of the SNR in the ED phase is largely de- 
termined by -B51, Mej and the power law index n of the 
ejecta envelope. The radius of the blast wave {Rb) during 
the ED phase is described by (from Eq. A5 of LH03) 



Rb 

Xo 



{n-2) 
ED 



(/)ED 47r(n — 3)n 

^ N (n-3)/(n-2) 



10 (n - 5) 
y (n - 3) 



tr 



(n-3)/2^ l/("-2) 



(1) 



where Zed = 1.19-f 8/^2 and 0ed = 0.39-0.6exp(-n/4) 
(the physical meanings of /ed and (/)ed are discussed in 
LH03 and TM99). Assuming E^i = 1, i?^, = 7.7 pc, nn 
= 0.1 cm~^, and t = 3000 yrs, the above solution gives a 
large ejecta mass of M^j ^ 20 Mq for the the power law 
index of the ejecta envelope n between 7 and 12. Adopt- 
ing a higher wind density of nn = 0.2 cm""^ reduces the 
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ejecta mass estimate and gives M^j 10 — 15 Mq. The 
ejecta mass is also reduced if we assume a less energetic 
explosion. With E^i = 0.5 and tt-h = 0.1 cm~'^, we es- 
timate Mej ^5 — 8 M0 . These estimates are based on 
the assumption that the remnant is strictly in the ED 
phase. The deceleration of the remnant during its tran- 
sition to the ST phase would increase the explosion en- 
ergy and decrease the ejecta mass. We note that, due to 
the high wind density near the center, radiative cooling 
may have been effective in the very early age of the rem- 
nant (which ma y be responsible for the optical emission 
from the eiecta. iLaming fc Hwand 1200 3*) and this could 
have affected the overall evolution. However, at this early 
time, most of the explosion energy would be in the form 
of the kinetic energy of the expanding ejecta, and the 
loss of energy by radiative cooling should be negligible 
in the total energy conservation and the overall remnant 
evolution. 

The shock velocity of the blast wave is given as 

V, = 2500 km s-^ ( ^] ( ( :-^)" 
\n-2j \7.7pcJ \3000yTsJ 

. For the observed electron temperature range of 0.6 
- 0.9 keV, this shock velocity indicates Te/Tp ~ 0.1, 
consistent with thos e in other remnants wit h compara- 
ble shock velocities ([Ghavamian "eFaI|[200l . The opti- 
cal observations of [O III] ejecta show an average expan- 
sion velocity of - ITOOkms"^ jGh avamian eFall [20051) . 
This velocity may be considered an upper limit for the 
current reverse shock velocity (the velocity in the rest 
frame of the ambient gas; ^). Using Eq. A13 of LM03, 
we derive ~ |^ near t = to. The derived velocity 

is insensitive to t because of its logarithmic depen- 
dency. We estimate ~ 2040 km s^^ for E51 1 and 

~ 1660 km s^^ for E^i ~ 0.5, which are consistent 
with the optical observations (jGhavamian "eFall [2OO5I) . 

While we have assumed a spherically symmet- 
ric explosion, recent observations and hydrodynamic 
modelings of SNe revealed a compelling case that 
high-ma ss SN explosions are i ntrinsically aspherical 
events (iWang fc Wheeled [200I iBurrows et all [20051: 
iJanka et al.ll2005[ ). Previous observations of G292.0-hl.8 
indeed suggested that t he explosion , of G2 92.0-I-1.8 may 
have been asymmetric. iPark et all (j2007| ) showed that 
the ejecta temperature is significantly higher generally in 
the northwestern and western regions of the SNR than 
in the southeastern regions. They suggested that this 
large-scale nonuniform distribution of the ejecta temper- 
ature may be caused by an asymmetric SN explosion: the 
explosion could have been more energetic toward NW-W 
than in SE. Winkler ct al. (2009) found that [O III] knots 
in north-south (N-S) axis are faster than those along the 
east-west axis, and suggested a possibility that explo- 
sions were more energetic along N-S axis. Also, there is 
a possibility that the azimuthal density variation around 
the remnant might be related with the asymmetric ex- 
plosion (§[31). However, these scenarios are not fully con- 
sistent with each other and further studies are needed. 
Our future study, which will focus on the X-ray charac- 
teristics of the shocked ejecta material, will help us reveal 
the detailed nature of the explosion of G292.0-I-1.8. 



5.3. Progenitor Star 

For the progenitor mass of ~ 20 — 40 M( t, estimated 
from previous X-ray observations (iHu ghcs & S inghl[T99^ 
[Gonzalez fc Safi-Harb|[200l iPark et al . 2004), the over- 
all observed wind properties are in plausible agreement 
with the general physical picture that the progenitor was 
a massive star which lost a significant amount of the mass 
through the RSG winds. However, given the uncertain- 
ties in the observational results and theoretical models, 
the current results do not provide a unique evolutionary 
track for the progenitor star (of a specific initial mass) 
which led to the explosion of G292.0-I-1.8. Based on our 
results, we briefly discuss a few possible scenarios for the 
evolutionary track of the progenitor. 

The progenitor star of initial mass of as low as ^ 20 Mq 
seems to be consistent with lower end of our wind mass 
estimate, where the star had a mass loss of ^ 10 — 15 Mq 
as the RSG wind while the rest of the mass went into 
the ejecta and the neutron star. In this scenario, the 
kinematics of the remnant can be explained by an explo- 
sion of i?5i ~ 0.5. While the stellar evolutionary models 
for a single star of M ~ 20 seem to predict smaller 
mass losses than 10 Mq, a possible companion star may 
have incre ased the mass loss , as proposed for the case 
of Gas A ([Young et al.|[20d6h . A significantly increased 
amount of mass loss is expected if the progenitor is more 
massive. For example, the models suggest that a star 
of 30 M0 loses ^ 20 M0 as the RSG wind ari d than ex- 
plodes as the RSG star ([Wooslev et al.|[2002f) . While a 
black hole might be expected after the explosion of such 
a massive star, the final fate of the compact remnant is 
inconclusive as observations show that a neutron star can 
be still formed fro m even more massive stars (> 40 Mq, 
[Muno et an[2006[ ). A progenitor star of > 35 Mq likely 
has evolved into the WR phase after the RSG phase. 
This might be the case of G292.0-f 1.8 if its WR phase 
has been short. In summary, given that the wind is most 
likely the RSG wind, the evolutionary models seems to 
prefer a progenitor's initial mass between ~ 20 — 35 Mq. 
This is plausibly consistent with the lower range of our 
wind mass estimates 15 Mq), although the higher 
mass progenitor cannot be completely ruled out. 

6. SUMMARY 

Using our deep ('^500 ks) Chandra observations, we 
present the results from our detailed analysis of the 
blast wave in the Galactic 0-rich SNR G292.0-H1.8. We 
find that the observed emission measure of the shocked 
ambient gas increases toward the contact discontinuity 
while the observed gas temperature decreases. Com- 
parisons o£_tlieseobservational results to self-similar so- 
lutions ([Ghevalierl [1982f l reveal that the observation is 
best described by the SNR currently expanding within 
the medium with a radially decreasing density profile 
(p oc r^^), likely a wind from the massive progenitor star. 
The inferred ambient medium is well described by a slow 
RSG wind from the progenitor star with a mass loss rate 
of M = 2 — 5 X 10~^MQyr~^ assuming a wind veloc- 
ity of Vu, ~ lOkms"^. The ambient density is too high 
for Wolf-Rayet wind, whose fast wind velocity implies 
a mass loss rate order-of-magnitude higher than what 
is generally accepted. The total swept-up mass of the 
wind is estimated to be ~ 15 — 40 Mq, where the quoted 
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range is due to the observed azimuthal density variation. 
The wind mass estimate may also be affected by other 
systematic uncertainties such as the distance and geome- 
try. Because the ambient medium is likely the RSG wind 
from a progenitor of M ^ 20 — 35 Mq , the higher range 
of our wind mass estimates (~ 40 M©) may not be fa- 
vored by standard stellar evolution models, although it 
is not completely ruled out. The overall kinematics of 
G292.0-I-1.8 is plausibly described by the models for a 
remnant expanding inside the RSG wind. 

Our results provide direct observational evidence for 
the blast wave interacting with the CSM winds in 
G292.0+1.8, which is consistent with the core-collapse 
origin of this SNR as revealed based on the central pul- 
sar and the 0-rich nature of the metal-rich ejecta. While 
the density of the CSM wind is one of the fundamen- 
tal properties to study the nature of young core-collapse 



SNRs and their progenitor's late-stage evolution, a di- 
rect observational measurement has been difficult. Deep 
X-ray observations of young core-collapse SNRs provide 
an excellent opportunity for such a study. Our future 
studies of the nature of the SN ejecta in G292.0-I-1.8 will 
be helpful to improve our understanding of the late-stage 
evolution of massive stars and their subsequent explosion 
as supernovae. 
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